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[57] ABSTRACT 

A symbol lock detector for an incoming coherent digi- 
tal communication signal which utilizes a subcarrier 
modulated with binary symbol data, dk, and known 
symbol interval T by integrating binary values of the 
signal over nonoverlapping first and second intervals 
selected to be T/2, delaying the first integral an interval 
T/2, and either summing or multiplying the second 
integral with the first one that preceded it to form a 
value X*. That value is then averaged over a number M 
of symbol intervals to produce a static value Y. A sym- 
bol lock decision can then be made when the static 
value Y exceeds a threshold level 6. 

3 Claims, 6 Drawing Sheets 
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SYMBOL LOCK DETECTION IMPLEMENTED 

WITH NONOVERLAPPING INTEGRATION 
INTERVALS 

5 

ORIGIN OF INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract, and is sub- 
ject to the provisions of Public Law 96-517 (35 USC 
202) in which the contractor has elected not to retain 
title. 

This application is a continuation, application Ser. 

No. 07/934,078, filed Aug. 14, 1992, now abandoned. 

TECHNICAL FIELD . 15 

This invention relates to symbol lock detectors for 
coherent digital communication systems after carrier 
and subcarrier synchronization has been achieved in a 
phase lock loop, and more particularly to lock detectors 2Q 
which may be implemented with a single integrator for 
processing nonoverlapping symbol intervals of an in- 
coming signal. 

BACKGROUND ART 

25 

Deep space network (DSN) receivers currently 
under development use phase locked loops (PLLs) to 
track the carrier, subcarrier, and symbol phase. Like 
most coherent receivers, the DSN receivers rely on 
lock detectors to provide the symbol lock status of its 30 
PLLs. Since carrier, subcarrier, and symbol synchroni- 
zation need to be achieved before any meaningful sym- 
bol detection can be initiated, symbol lock detectors 
play a vital role in the final decision of accepting or 
rejecting the detected symbols. In the past, symbol lock 35 
detectors have employed overlapping symbol intervals 
in their operations and therefore require two integrators 
operating over staggered time intervals. 

During operation, a loop is assumed to be locked 
when its lock indicator consistently has a positive status. 40 
The carrier and subcarrier lock detectors currently used 
in DSN receivers have already been analyzed. The 
present invention concerns the analysis of three new 
symbol lock detectors which simplify implementation 
for DSN receivers and two prior-art symbol lock detec- 45 
tors for comparison. 

Statement of the Invention 

In accordance with the present invention, an incom- 
ing signal that has undergone carrier and subcarrier (if 50 
any) synchronization is processed to determine whether 
or not symbol lock has also been achieved. The symbol 
data, dfc takes on the binary values of a pulse p(t) of a 
known duration, such as in a nonretum-to-zero (NRZ) 
or Manchester code, where the probability of the data 55 
having a + 1 is equal to having a — 1 (or 0) value, (i.e., 
the data transition probability equals one half). While 
the prior-art lock detectors employ two integrators for 
processing two overlapping symbol intervals, the lock 
detectors of the present invention require only one 60 
integrator for processing two nonoverlapping time in- 
tervals, e.g., integration of the last half of one (d*.i/2) 
and the first half (d/c/2) of the next symbol interval, 
followed by integration of the second half (d*/2) of the 
next and the first half (dfc+ 1/2) of the following symbol 65 
interval, or the first half of each symbol interval fol- 
lowed by the second half of each symbol interval. The 
receiver is assumed to have perfect knowledge of the 
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symbol interval T=kt, where k is the number of data 
bits, and t is the data bit interval. 

The first integration output is processed and delayed 
a half symbol interval, T/2, for determining symbol 
lock by addition to or multiplication with the following 
integration output, where the processing of the first and 
second integration outputs is a squaring operation or an 
absolute value operation, and symbol lock is determined 
by addition of the processed second integration output 
to the processed first integration output delayed a half 
symbol interval and accumulating the sums. Symbol 
lock is then determined by repeatedly thresholding the 
accumulated sums every M symbol intervals. An alter- 
native arrangement for determining symbol lock is by 
multiplying the second integration outputs by the first 
integration output delayed a half symbol interval to 
obtain a signal power estimation, and after accumulat- 
ing signal power estimates over M symbol intervals, 
thresholding the accumulated estimates. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a functinal block diagram common to sym- 
bol lock detectors shown in FIGS. 2 a through 4. 

FIG. 2a is a functional block diagram of a prior-art 
square-law symbol lock detector with integration of 
overlapping intervals (SQOD). 

FIG. 2b is a functional block diagram of a prior-art 
absolute-value symbol lock detector with integration of 
overlapping intervals (AVOD). 

FIG. 3a is a functional block diagram of a square-law 
symbol lock detector with integration of nonoverlap- 
ping intervals (SQNOD) in accordance with the present 
invention. 

FIG. 3 b is a functional block diagram of an absolute- 
value symbol lock detector with integration of nono- 
verlapping intervals (AVNOD) in accordance with the 
present invention. 

FIG. 4 is a functional block diagram of a signal- 
power symbol lock estimator (SPED) with integration 
of nonoverlapping (first and second) halves of incoming 
signal symbols in accordance with the presnet inven- 
tion. 

FIG. 5a and 5b are graphs of the probability density 
function of the output Y for the SQOD of FIG. 2a when 
symbol lock is not present and the output Y has a high 
SNR =5 dB and a low SNR=— 5 dB, respectively. 

FIG. 6 is a graph of the probability of lock detection 
when t is an unknown constant over a decision interval 
for symbol lock detection of FIGS. 2a, 2b, 3a, 3b and 4. 

FIG. 7 is a graph of the probability of lock detection 
versus SNR when r is uniformly distributed and chang- 
ing from symbol to symbol for symbol lock detectors of 
FIGS. 2a, 2b, 3a, 3b and 4. 

FIG. 8 is a graph of the probability of lock detection 
versus SNR when the false alarm rate is computed in 
the absence of a signal for symbol lock detectors of 
FIGS. 2a, 2b, 3a, 3b and 4. 

DETAILED DESCRIPTION OF THE 
INVENTION 

The symbol lock detectors considered are divided 
into two groups. The detectors in the first group are 
prior-art symbol lock detectors that process the over- 
lapping outputs of two symbol data bit integrators, 
whereas those in the second group, which depart from 
the prior art, use one integrator for two nonoverlapping 
intervals of an incoming signal to determine whether or 
not symbol lock has been achieved. 
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The first group of symbol lock detectors employ 
either the squares of the two overlapping integrator 
outputs or the absolute value of the integrator outputs. 
The second group consists of three new lock detectors, 
two of which use the same mathematical operations of 5 
squaring or taking the absolute values of integration 
outputs of nonoverlapping half symbol intervals but 
only one integrator, while the third detector, which also 
uses only one integrator, does not include either squar- 
ing or absolute value operations, and instead functions 10 
as a signal power estimator by multiplying the two 
integration outputs of nonoverlapping half symbol in- 
tervals. 

The five lock detectors are compared based on the 
lock-detection probability as a function of the symbol 15 
SNR for a given false-alarm probability and a fixed 
observation interval. 

Although symbol lock detection has been addressed 
before [J. K. Holmes, Coherent Spread Spectrum System, 
New York: John Wiley and Sons, 1982 and K. T. Woo, 20 
Shuttle Bit Synch Lock Detector Performance, TRW IOC 
No. SCTE 50-76-1 84/KT W,TR W Corporation, ‘ El 
Segundo, Calif., Apr. 5, 1976], the analyses have ne- 
glected the interdependence between symbol synchro- 
nizer bandwidth, and lock detector bandwidth. The 25 
symbol synchronizer bandwith refers to the one-sided 
loop noise bandwidth Bl of the digital data- transition 
tracking loop [M. Simon, “An Analysis of the Steady- 
State Phase Noise Performance of a Digital Data-Tran- 
sition Tracking Loop,” JPL Space Programs Summary, 30 
37-55, Vol. 3, Jet Propulsion Laboratory, Pasadena, 
Calif., pp. 54-62, Feb. 28, 1969] used in the receivers. 
The lock detector bandwidth is defined as the fre- 
quency at which the lock detector provides a status 
signal of being in or out-of-lock. For example, the lock 35 
detectors considered here indicate loop status once 
every M symbols. Consequently, the bandwidth of 
these detectors is 1/MT, where T is the symbol interval. 

The probability of false alarm, P/ 0 , is defined in two 
ways. In the classical sense, it is defined as the probabil- 40 
ity of declaring a signal (or target as in radar applica- 
tions) to be present when it is not present. In deep space 
applications however, it is more appropriate to define 
P/a as the probability of declaring a loop to be in-lock 
when it is out-of-lock. That is, declaring the timing 45 
error to be zero (inlock) when the loop is slipping cycles 
and operating with a non-zero timing error (out-of- 
lock). 

In discussion below, the false alarm rate is shown to 
be drastically different depending on the definition 50 
used. In addition when the loop is slipping cycles, the 
false alarm rate is shown to depend strongly on the ratio 
of the lock detector bandwidth to the symbol loop 
bandwidth. For example, when the loop is slipping and 
1/B l < <MT the lock detectors operate with accept- 55 
able false alarm rates because there are several uncorre- 
lated samples of the timing error 1 /t within the MT 
second decision interval. On the other hand, when 
l/Bf.>>MT the false alarm rates are unacceptable 
because the timing error is constant over several deci- 60 
sion intervals. Note that a good rule of thumb is to 
assume that the loop provides uncorrelated phase esti- 
mates every 1/Bi seconds. As a result, the symbol tim- 
ing error at time t,- is uncorrelated with the symbol 
timing error at time t y - when [t/ — 1/[ = 1/Bi. This article 65 
considers the special cases of l/Bi=MT and l/Bi=T. 
The first case is analyzed and simulated whereas, the 
second is simulated but not analyzed. When the thresh- 
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old is adjusted in the presence of noise only, the perfor- 
mance can be derived from the previous analysis by 
setting the signal amplitude to zero. 

Generic Description of Lock Detectors 

FIG. 1 is a block diagram showing the signal process- 
ing functions common to the symbol lock detectors 
analyzed and discussed below, including both prior-art 
lock detectors of FIGS, la and 2b, and new lock detec- 
tors in accordance with the present invention as shown 
in FIGS. 3 a, 3b, and 4. The received signal is assumed to 
have been mixed with perfect carrier and subcarrier 
local reference signals so that the input to the lock de- 
tectors is a baseband signal of the form 

dt)=Ad(t)+n(t) (1) 

where A is the signal amplitude and A 2 is the received 
data power with 

00 ( 2 ) 

4(f) =2 d^f-kT) 

*= — oo 

and where n(t) is the additive white Gaussian noise 
process with a single-sided power spectral density 
(PSD) No (Watts/Hz). The data symbol d* takes on the 
binary value -f 1 and — 1 (or 0) with equal probability, 
and p(t) is the received data pulse shape of duration T 
seconds. For comparison purposes, only the NRZ pulse 
is considered in the analysis but the results can be ex- 
tended to any pulse shape, such as Manchester encoded 
data. The receiver is assumed to have perfect knowl- 
edge of T, but not the symbol epoch, i.e., the receiver 
has estimated perfectly the symbol rate but not neces- 
sarily the start and end of the symbols. 

The signal processing functions for the symbol lock 
detector 10 in FIG. i depend on the processing of inte- 
gration outputs in the detector. Its output X* is at the 
symbol rate and typically many samples of X* are aver- 
aged in an accumulator 12 to obtain the decision statistic 
Y compared in a block 14 with a threshold value 8. If Y 
is greater than the threshold 8, the loop is declared to be 
in-lock, otherwise it is declared to be out-of-lock. 

The symbol timing error (parameter) r in FIG. 1 is 
the phase error between the symbol phase and the phase 
estimate provided by the symbol synchronizer. The 
in-lock case is analyzed by setting the timing error r to 
zero. In practice, the error is not identically zero, but it 
is a very small value. When there is a signal present, the 
out-of-lock model for r depends on the relation between 
Bl and 1/(MT). When Bl= 1/(MT), t is modeled as an 
unknown constant over a decision interval (MT sec- 
onds) but independent and uniformly distributed from 
one decision interval to the next. 

Alternatively, when Bl= IT, the timing error is mod- 
eled as constant over a symbol interval, T, but indepen- 
dent and uniformly distributed from symbol to symbol. 
In this case, if the decision time MT> >T (as it usually 
is), each decision statistic encompasses the entire range 
of t. When there is no signal present, the model of r is 
irrelevant because the out-of-lock performance is inde- 
pendent of T. 

The respective prior-art and new detectors consid- 
ered and compared below are the square-law detector 
with overlapping (SQOD) and non-overlapping 
(SQNOD) integrators shown in FIGS. 2a and 3a, the 
absolute-value detectors with overlapping (AVOD) 
and non-overlapping (AVNOD) integrators shown in 
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FIGS. 2b and 3b, and finally the new signal power 
estimator detector (SPED) shown in FIG. 4. 

In the prior-art SQOD detector, FIG. 2a, the input 
signal r(t) is integrated over two symbol periods by two 
integrators 21 and 22: one in phase with the estimated 5 
symbol interval and the other staggered by half a sym- 
bol duration. The resulting inphase and quadraphase 
samples I* and Q/tare correlated due to the overlapping 
intervals. The quadrature samples I* are squared in a 
processor 23 and. delayed a half symbol period by a 10 
delay element 24 while the samples Q k are squared in a 
processor 25. The delayed (I*) 2 and the undelayed (Qk) 2 
samples are then combined in a summing circuit 26 to 
form an output X* which are averaged in an accumula- 
tor 27 to obtain the decision statistic Y which is then 15 
compared with a predetermined threshold 8 in a detec- 
tor 28 to reach a lock decision. The prior-art AVOD 
detector, FIG. 2b merely replaces the squaring opera- 
tions in FIG. 2a with absolute value operations in 
blocks 23' and 25'. Hence, the I and Q samples of the 20 
prior-art AVOD detector are also correlated. 

The new SQNOD detector, FIG. 3a, processes the 
integration outputs of nonoverlapping half symbol in- 
tervals from a single integrator 30. A single squaring 
operation in block 31 produces the square of the inte- 25 
gration outputs from both halves of a symbol interval, 
and a multiplexer 32 separates the first and second half 
symbol interval integration outputs sending the first half 
through a delay element of half a symbol interval and 
the second half directly to a summing circuit 34. Its 30 
output Xk is then averaged in an accumulator 35 and its 
output, and the decision statistic Y is compared with a 
threshold 8 in a detector 36. As before, replacing the 
squaring operations in FIG. 3 a with absolute value 
operations in block 31' yields its counterpart the 35 
AVNOD detector shown in FIG. 3b. The integrator 
outputs in these cases are uncorrelated because the inte- 
grated intervals are nonoverlapping. 

The SPED detector shown in FIG. 4, which also uses 
a single integrator 40, was considered for symbol lock 40 
detection because it already existed as part of a split 
symbol moment SNR estimator [K. T. Woo, supra] in 
DSN receivers used for SNR estimation. The inphase 
(I*) and quadraphase (Qk) integration outputs of the 
SPED are obtained by integrating the received signal 45 
over the respective first and second halves of a symbol 
interval and separating them by a multiplexer 41. Since 
the noise in the first and second half are independent, 
delaying the inphase integration output, lk, in a delay 
element 42 and then forming the product of I* and Qk in 50 
a multiplier 43 and averaging over number, M, of sym- 
bol intervals in an accumulator 44, provides as a deci- 
sion statistic Y an estimate proportional to signal power 
which is then compared to the value 8 in a threshold 
detector 45 to provide a lock decision. 55 

Performance Analysis 

In-lock performance is measured in terms of the prob- 
ability of declaring a phase-locked loop (PLL) as being 
symbol locked when there is no timing error, i.e., the 60 
probability that the decision statistic Y is greater than 
the threshold value 8 when t= 0. Note that r=0, or no 
phase tracking error, is equivalent to setting the symbol 
synchronizer loop SNR to infinity. The degradation in 
detection probability due to timing jitter (non-infinite 65 
loop SNR) is minimal and has been addressed in the 
case of carrier lock detectors [A. Mileant and S. Hinedi, 
“Costas Loop Lock Detection in the Advanced Re- 
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ceiver,” TDA Progress Report 42-99, Vol. July-Sep- 
tember 1989, Jet Propulsion Laboratory, Pasadena, 
Calif, pp. 72-89, Nov. 15, 1989]. The out-of-lock perfor- 
mance is measured by the probability of false alarm, i.e., 
the probability of declaring the loop as locked when it 
is not locked. The out-of-lock performance in the pres- 
ence of a signal is analyzed for the case Bi=l/(MT). 
Note that, in that case, the timing error r is independent 
from one symbol to another and the decision is per- 
formed after averaging a number, M, of symbols. On the 
other hand when B^= 1/(MT), r is an unknown con- 
stant during a decision interval and independent from 
one decision to the next. Setting M = 1 in the latter case 
would imply a decision every symbol, which is funda- 
mentally different from the case Bz.= 1/T decision 
using M symbols. Hence, the performance when 
B£= 1/T cannot be derived from the case of 
Bl= 1/(MT), simply by setting M=l. The out-of-lock 
performance when there is no signal present is also 
analyzed. 

In this discussion, only the final equations are shown. 
Derivations of these equations are set forth in the vari- 
ous appendices A through F for all five symbol lock 
detectors SQOD, SQNOD, AVOD, AVNOD and 
SPED which, by this reference, are hereby made a part 
hereof. In all cases, the decision statistic can be ex- 
pressed as the average of samples Xk over M symbol 
intervals given by the equation 


l 

M 


M 

■ 2 X k 
k= 1 


(3) 


Note that the random variable X* is peculiar to each 
detector. When the timing offset r=0, the adjacent 
samples Xk and Xk+ 1 are correlated in the two prior-art 
detectors (SQOD and SQNOD), whereas for the three 
new detectors (AVOD, AVNOD and SPED) they are 
uncorrelated. In all cases, the random variable X*is not 
Gaussian due the nonlinear operations on I* and Q*. For 
large M, the random variable Y is modeled as Gaussian 
due to the Central Limit Theorem (CLT). The CLT 
theorem applies to the sum of correlated random vari- 
ables when none of the variables being summed or mul- 
tiplied dominates over the others [D. Fraser, Non-Para- 
metric Methods in Statistics, New York: John Wiley and 
Sons, 1957]. This model for Y is justified by simulation 
results. The probability of lock detection is the probabil- 
ity that the Gaussian random variable Y surpasses the 
threshold 8. Hence, it is given by 


Pd 


-J 


8 N 2 7r<T}2 


exp 


r (y - PY) 2 1 , 

L — 


<“> 


where fry and ay 2 are the mean and variance of Y when 
r is exactly zero. Using the definition of the error func- 
tion 


erflx) = -rp- f exp 

Ntt J 0 


(5) 


one has 


Pd = \ Y erf 




MT 


cry 


( 6 ) 
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1 1 , 8 

~-~ erf {^r, 



r (k + 1 )T + T 

h = Kt)dt 

J kT+ t 

= d k A(T - r) + + iVK*) + Ni(k) 

and the quadraphase integrator output is given by 


where SNR# denotes the detector SNR defined by 10 n _ f ^ + yx>T + T , 


■‘te-)’ 


f VA. T J/A.JX T T 

Qk = Hfidt 

J (k + i)T + r 


dkA dk+\A + T 


The threshold S is chosen to maintain a fixed proba- 
bility of false alarm. The probability of false alarm is the 
probability that the out-of-lock decision statistics do not 
surpass the threshold. Hence, it is given by 20 


N 2 (k) + Ni(k + 1) 


0St< y 


dk+\A [~Y~ - T J+ dk+iA - y-^H 


N 2 (k) + N\(k + 1) 




where r is limited to the interval [O.T], and where 


where /<>( Y) is the out-of-lock density of Y. The thresh- 
old 8 is computed by solving Eq. (9) for a fixed P/ a . 
When there is a signal present and B/.= 1/(MT), the 
statistic Y is no longer Gaussian and f 0 ( Y) must be 30 
obtained numerically or by simulation as shown below. 
When there is no signal present, the CLT can be in- 
voked and the out-of-lock decision statistic can be mod- 
eled as Gaussian. This model is verified by simulations 

-JC 

below. In this case, Eq. (9) can be written as 


1 1 J 8 - y-Yo 

p f“- 2 ~ 2 er f \rr 

^ N2 c t Yo 


where fiYoand 0- re the out-of-lock mean and variance of 
the decision statistic Y. The threshold 8 is given by 


r (k + i)T + r 

Ni(k) = n(t)dt 

J kT 4- r 


Ak+\)T+T (16) 

N 2 (K) = n(i)dt 

J (k + i)T+r 

Since n(t) is a white Gaussian process with one-sided 
PSD No, the N,’s are independent Gaussian random 
variables with mean zero and variance cr„ 2 =(N 0 T)/4. 
The samples X*=Lfc 2 — Q* 2 and summing M of them 
yields Y. From Appendix A, the in-lock mean and vari- 
ance of Y are given by 


VsNoT 
MY = J 


°ToT + HYo 


where Y=erf- 1 (l — 2Pf a ). Substituting Eq. (11) into 
Eq. (6) relates the probability of detection to the no 
signal (classical) false alarm rate, namely, 




VT <TYo Y + I ‘■Yo - M-X 

VT7 r 


r ^=(^¥-) _^(-x-+ 2 %+t)- 


2 (^-i)(-T-t)_ 


The out-of-lock mean and variance when there is a 
signal present and r is an unknown constant over a 
decision interval (B#=1/(MT)) are given by 


The next five subsections derive the in-lock and out- 
of-lock mean and variance for all five detectors SQOD, 60 
SQNOD, AVOD, AVNOD and SPED for compari- 
son. 

(1). Square-Law Lock Detector with Overlapping Inter 
vals (SQOD) 65 

The SQOD detector is shown in FIG. 2a. For the 
input given by Eq. (1), the inphase integrator output is 
given by 


(■ nsir - tt - t )■ * <" - «" - » -rr ] 


where % denotes symbol signal-to-noise ratio and is 
defined as 
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% = 


( 21 ) 


10 

-continued 


No 


Setting A=0 in Eq. (20) and substituting the result in 
Eqs. (19) and (20) yields the out-of-lock mean and vari- 
ance in the no signal case. Hence, the no signal mean is 
zero but the no signal variance is given by 


yar^X/c) = (NoT)^ ~Y~ + 1 — Vs Il(Vs> + Gl(%) + 


(29) 


2/l(i?,) + 2/ 2 (t),) + m-n: 


I,) ^ 


10 


o = 0) = (^) [m- f - 2(M -1)4.] } Co^X k JC k+i ) = (Ar 0 7){ - e - x - p( 2g 2%) + (f + 4-)> 


(30) 


(2). Absolute- Value Lock Detector with Overlapping 15 / 

Intervals (AVOD) er/(NV , + exp (-^W ) 

For the AVOD detector shown in FIG. 2b, absolute ^ ' 

values are used instead of squares. The expressions for 
I* and Qyt given by Eqs. (13) and (14) are still valid but 
now X*= 1 1*1 — |Q* I - From Appendix B, the in-lock 20 % 

mean and variance are given by 


8 Not), ) 

J- I 2=3(7),) + ( 73 ( 7 ),) + 2 / 4 ( 7 ),) - G 4 (Vs) - 2 / 5 ( 7 ?,) - 


• 2 




=WT) 


exp(-T),) - 1 


N7 


'+’11?, 


»/(N 7), 


(23) 


25 


and 


o-y 2 = -4j- [M VariXk) + M(M - \)Cov(X k Xk + \)} 


(24) 


Coy 0 (X k ,X k+ j) 


Y 2=*i(t),) - Gi(t?,) - 2G2(7),) - 22/2(7),) ] 'j 

,(«*) 


30 


where 


I [2=3(1?,) + G 3 (t),) + 2/4(7),) ■ 


204(1,,) - 22/5(7),)] for/ is 2 


^/O*) = (2V< 


r o7>{' 


3% 


1 - yAF\(Vs) + 2=2(7,,)] 


(25) The functions F/, G, and H, in Eqs. (28) through (31) 
are defined in Appendix B and plotted in FIG. B.l 


[exp(— 2%) - 2exp(— 7,,) + 1] - ^f- erfy N t), ■ 


4ir 


35 


40 


versus %. Setting %=0 in Eqs. (28) through (31) yields 
the out-of-lock statistics in the no signal case. The no 
signal mean is zero but the variance is given by 


sfT ,*42 _ a ) ^ . (A£)[„ + J( „ _ „ ( J_)] 


(32) 


and 


Cov(X k J k+ i) = (N c 


r 0 7){ 


exp(-2i),) + exp(-i),) 
2n 


(,g) 45 (3)- Square-Law Lock Detector with Nonoverlapping 
Intervals (SQNOD) 

The SQNOD detector is shown in FIG. 2Z?.For the 
input of Eq. (1), The inphase and quadraphase integra- 


Vs_ 

2 


erft'fvT ) + [exp(— Tjj) + j-J^NT^T ) X 50 t0r OUt P UtS are gi ven b y 


N 7 T 2 


[Fi( % ) + F 2 (Vs)J 


) 


ik 


s 


(* + 2 )T+t 
(k + i)T + r 


*t)dt 


(33) 


55 


The out-of-lock mean and variance when there is a 
signal and r0 is constant over M symbols are given by 


V-Yo = 0-° 
and 

_2 _J 

M 2 

where 


(27) 


60 


dk-^Y~ + 2Vl(*) 


d k A 


T) 


(28) 


(Af - 1XA2 - 2)CovJiX k Jk+2)] 


65 


and 


Qk = j 


d k + 1 A I t 


(k + \)T+r 
(k - i)T + t 


0 = T < " 


T T 

+ 2V,(*) 


<l)dt 


(34) 
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ry 2 = 


12 

-continued 


( N 0 T 

^ M 


^Hr + l- -iir[ I + 5ex p ( -^ 


(41) 


11 

-continued 
dk-\A - T^j + 

~ dkA^Jf- + T j+A'iW 0^r<-y- 
d k -¥- + N 2 (k) L. St< Z. 

The noises Ni(k) and N 2 (k) are given by Eqs. (15) and 
(16) after changing the integration limits to those in Eqs. 
(33) and (34) . As a result, they are independent Gauss- 
ian random variables with zero mean and variance 
cr n 2 =(N 0 T)/4. The sample X* is the difference of the 
squares (i.e. , X*=I* 2 — Q* 2 ). From Appendix D, the 
in-lock mean and variance of Y are given by 

flsNoT (35) 

nr = — T~ 


[ 5ex P (— t)+ 1 ]- nr er/ j 

15 The out-of-lock mean and variance when there is a 
signal present and r an unknown constant over MT 
seconds is given by 

2Q pr 0 = 0.0 (42) 

and 


and 




-(4F-X 


i ), 2 + 24% + 16 
64 


(36) 25 


For the case of false lock with signal present and r an 
unknown constant over M symbols, one obtains 30 


°- 2 Yo 


PYo = 0.0 


(37) 

and 



f Nj-T 1 \ 


(38) 

i! 

— 1 S "x 

( % 2 + 25% + 15 

) + M(M-l)-^L] 


\ 60 


35 


40 


cPr 0 = -jp [M VarJX k ) + M(M - \)Cov 0 (X kl X K+i )] 


(43) 


Var 0 (X k ) = (Nol) 



47T 


• exp(-%) ■ 


(44) 



When there is no signal present the out-of-lock mean is 
zero and the variance is given by setting %=0 in Eq. 
(38) . Consequently, the no signal out-of-lock variance 45 
is 



•to-.-(4£ )0) 

(4). Absolute- Value Lock Detector with Nonoverlap- 
ping Intervals (AVNOD) 

This detector shown in FIG. 3 b is the same as the 
SQNOD detector shown in FIG. 3 a with the squaring 
operations replaced by absolute value operations. 
Hence, Eqs. (32) and (33) for I u. and Q& are valid but 
now Xk= 1 1* | — | Q* | . From Appendix E, the in-lock 
statistics for Y are given by 


(40) 



where the function Z is defined in Appendix E and 
plotted in FIG. B.2. 

For the out-of-lock case with no signal, the mean is 
50 zero and the variance is obtained by setting 7] s =0 in 
Eqs. (43) through (45). Hence, the out-of-lock variance 
is given by 


55 


YoiVs = 0) = 


N 0 T 

M 


^ 2 4 it ) 


(46) 


(5). Signal-Power Estimator Lock Detector (SPED) 

In the detector shown in FIG. 4 , denote the integra- 
tions over the first half of the assumed symbol interval 
as I* and the second half as Qh Then, the I* and Q* 
samples are given by 


65 


Ik = d k ~~ + mk) 


(47) 


and 


and 
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-continued 


: = rf *^( T-r) 


+ + AWO 


and X*=I*Qi. From Appendix F, the in-lock mean and 
variance of Y are 




The out-of-lock case with signal present has mean and 
variance when r is constant over M symbols is given by 


, X ..( 5% 2 + 20 % + 12 ^ % 2 "I 

<r ™ M 2 l M \ 192 J+M(Af 1) 192 J 

As before, the out-of-lock variance in the no signal case 
is given by 


Discussion and Simulation Results 

Digital simulation was used to verify the foregoing 40 
analysis. In the out-of-lock state for a long-time constant 
Bi=l/(MT), the symbol timing error r is modeled as 
constant over a decision interval (MT seconds) but 
independent and uniformly distributed over the collec- 
tion of all decision intervals. The timing error in the 45 
in-lock state is modeled as being zero. Although the 
special case of r constant over M symbols was analyzed 
for performance comparison purposes, it is not advis- 
able to operate a practical system under these condi- 
tions due to unacceptable false alarm rates. This case 50 
has higher than usual false alarm rates because the deci- 
sion statistic for small values of r is not significantly 
different from the statistic for t=0. As a result, the 
out-of-lock states corresponding to small values of r are 
frequently declared to be in-lock because they are mis- 55 
taken for the case when r=0. This problem can be 
ameliorated by lengthening the observation time rela- 
tive to the time constant of r (i.e., shortening the time 
constant of t). In practice, it is recommended that the 
observation time be at least ten times longer than the 60 
time constant of t. As noted above, the out-of-lock 
density function for Y in this case is not Gaussian. Con- 
sider the decision statistic Y when the loop is out of lock 
and r is constant over M symbols. In general, it can be 
written as 65 

Xk=siAr)+nk+SkiT)nk (54) 


where, in all five detectors, the signal term s* is random 
and uniformly distributed because r is a uniformly dis- 
tributed random variable. The density of the noise n* 
depends on the detector being implemented. Summing 
M samples of X& where Xk is at the symbol rate in all 
cases, yields the decision variable Y. Since r is constant 
over the sum, at high SNR (i.e., for strong signal levels) 
the density function of Y approaches a uniform distribu- 
tion as shown in FIG. 5a. However, at low SNR the 
noise term dominates and the density of Y is Gaussian 
due to the central limit theorem as shown in FIG. 5 b. 
The density in FIGS. 5 a and 5 b was obtained via nu- 
merical integration as well as simulation. Both methods 
are seen to agree very well. The numerical method 
computed the density function of Y, /( y), by averaging 
over t the conditional pdf /( y/r). The latter is Gaussian 
with mean and variance both a function of r. The simu- 
lation method computed the histogram of Y and then 
set f(y) — [P(y — Ag Yi£y+ A)]/A, where A is the size of 
a histogram bin. The histograms were generated using 
1,00,000 million symbols which corresponds to 10,000 
decisions (Y’s), since there are 100 symbols/decision. 

FIG. 6 compares the probability of detection perfor- 
mance of all five detectors for M=100 and P/ O =0.25. 
Note that the overlapping detectors SQOD and 
AVOD, which are identical except for the squaring and 
absolute value operations, have nearly identical perfor- 
mance. As expected, the AVOD is slightly better at 
high SNR, whereas the SQOD is slightly better at low 
SNR. The non-overlapping detectors SQNOD and 
AVNOD also have nearly equal performance. Once 
again, the absolute value operation yield better results at 
higher SNRs. The signal level estimator (SPED) is 
better than the non-overlapping detectors but worse 
than the overlapping detectors. The probability of de- 
tection results in FIG. 6 change when P f a or M change. 
For example, increasing the observation interval in- 
creases the detection probability because it increases the 
detector SNR = (p- y 2 / cr yu 2 ). Accepting a higher false 
alarm rate increases the probability of detection because 
it lower the threshold 8. In generating these curves, 
50, (XX) symbols were simulated for each value of SNR. 
Since there are 100 symbols/decision, the detection 
probability for a given SNR is based on 500 decisions. 

In the out-of-lock state for a short-time constant 
(Bl= 1/T), the symbol timing error z is modeled as 
uniformly distributed random variable that changes 
independently from symbol-to-symbol. For this case, 
the probabilities of detection for all five detectors are 
computed by simulation for M= 100, Pf a = 10 _2 and the 
threshold S is set according to Eq. (1 1). The false alarm 
rate was verified by simulation. The results are plotted 
in FIG. 7 versus symbol energy-to-noise ratio 7 ] s . In 
these computer simulations, the detection probability 
for a given SNR is based on 40,000 decisions. 

The results show that, the AVOD performs slightly 
better than SQOD at high SNR, whereas they seem to 
perform identically at low SNR. The nonoverlapping 
detectors SQNOD and AVNOD also have nearly equal 
performance at low SNR, but AVNOD performs about 
1 dB better for values of symbol SNR higher than —4 
dB. As far as the SPED, it performs about 2 dB worse 
than the overlapping detectors and 3 dB better than the 
other two nonoverlapping detectors. Also by simula- 
tion, the false-alarm rate that was used in setting the 
threshold was verified. 

In the situation of no signal, i.e., when there is no 
signal present as distinguished from the case when there 
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is a signal and r=0, the out-of-lock statistic is Gaussian 
with zero mean and the in-lock statistic is Gaussian with 
non-zero mean. Probability of detection results are com- 
pared in FIG. 8. Interestingly, the performance of the 
overlapping and nonoverlapping detectors are grouped 5 
together, but the signal level detector (SPED) now has 
the best performance. The interdependence between P^, 

P f a and M is the same as in the other two cases. 

Conclusion ]0 

The performance of two prior-art symbol lock detec- 
tors shown in FIGS, la and lb have been compared 
with three new ones shown in FIGS. 3 a, 3b and 4 . They 
are the square-law detector with overlapping (SQOD) 
compared with nonoverlapping (SQNOD) integrators 15 
and others, the absolute value detectors with overlap- 
ping (AVOD) compared with nonoverlapping (AV- 
NOD) integrators and others, and the signal power 
estimator detector (SPED) compared with all others. 
The analysis considered various scenarios when the 20 
observation interval is much larger or equal to the sym- 
bol synchronizer loop bandwidth, which has not been 
considered in previous analyses. Also, the case of 
threshold setting in the absence of signal was consid- 
ered. The analysis has shown that the square-law detec- 25 
tor with overlapping integrators (SQOD) outperforms 
all others when the threshold is set in the presence of a 
signal, independent of the relationship between loop 
bandwidth and observation period. The square-law 
detector and absolute-value detector with overlapping 30 
integrators outperformed corresponding detectors with 
nonoverlapping integrators, but implementation of the 
SQOD and AVOD require two separate integrators 
since both integrators must operate at the same time due 
to overlapping, whereas the SQNOD and AVNOD 35 
symbol lock detectors may use a single interval since 
the separate integration operations required do not 
overlap in time. On the other hand, the signal-power 
estimator detector (SPED) outperforms all others when 
the threshold 8 is set in the presence of noise only, and 40 
it requires only a single integrator for implementation. 

Appendix A 

Derivation of the Mean and Variance of the SQOD 

The inphase and quadraphase integrator outputs are 45 
given by Eqs. (13) and (14), respectively. The output of 
the lock detector X,t=I* 2 — Q* 2 . Consequently, 


p.k-e{Ik 2 }-<i{Qk 2 } 
where p.k=e{Xjt}. 


Vc^Xk) = <={[* 2 - fit 2 ] 2 } - e 2 {/* 2 - Q k 2 } 

= C W + Qk A - 2Ik 2 Qk 2 > - M 2 

The covariance of X*, with Xk+j is 

CovQLhXk+j = c\l ] “1 -f~ tiQlP'Qk+fl ~ 
t{Ik 2 Qk+j 2 } - dlk+j 2 Qk 2 }-ILkHk+j (A-3) 

When the loop is in lock, Eqs. (A-l) through (A-3) are 
evaluated with t set to zero in Eqs. (A- 10) through 
(A-18). Hence, the in-lock moments of X* are given by 


Var{X k ) = ■ 


+ 8A 2 T 2 cr n 2 + 12 <r„ 2 


Cov(X k Jk+ 1 ) = -2A 2 T 2 cr n 2 - 2cr„ 4 (A-5) 

and for j=2, this can be shown to be 

Cov(X k ,X k+ j )= 0 (A-6) 

When the loop is out-of-lock, r is modelled as a uni- 
form random variable. Using this model for r in Eqs. 
(A- 10) through (A-18) and substituting the results into 
Eqs. (A-l) through (A-3) give the out-of-lock moments 
of Xk . Namely, (where the additional subscript 0 de- 
notes out-of-lock), 

V aro (Ik 2 - Qk 2 ) = pyp + A 2 T 2 cr 2 + 120 -/ ( } 

Cm 0 (XkJC k+ 1) = - A 2 T 2 a 2 - 2o-/ ( ’ 

and for j =2, this can be shown to be 


Cov 0 (X k Xk+j) = pp 


The following equations were used to compute the 
variance of X k and covariance of X* with Xfc+i: 


t 2 } = A 2 (T 2 — 27€{t} + 2«{t 2 }) + 2 o-„ 2 


<Kie* 2 |} = 


A 2 (—. : — eifl} + 2 £i(t 2 }) + 2o-„ 2 ei{l} 


0 = T < ■ 


A 2 (-2j— £ 2 {1} - 47lE2{r} + 2£2{t 2 }) + 2o-„ 2 €2{1) y- = t < T 


e{/* 4 } = ^ 4 (r* - 47 Mr} + 12r 2 <d> 2 } - 167e{T 3 } + 8el> 4 }) + ^ 2 o-„ 2 (12r 2 - 24Te{r} + 24<={r 2 }) + 12o-„ 4 


cie/} = 


A\-±- £1 {1} + 8 £ ,{t 4 }) + 

A 2 o-„ 2 (6T 2 e]{\} + 24ei{r 2 }) + 12 ct„ 4 £ i{1} 
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-continued 

A\ -^~ - £ 2 {1} - 327 3 £ 2 {t} + 487- 2 c 2 {t 2 } - 327-£ 2 {t 3 } + 8£ 2 {t 4 }) + 

/f 2 or„ 2 (30r 2 e 2 {l} - 487£ 2 {r} + 24£ 2 {t 2 }) + 12<r„ 4 e 2 {l} -j- = r < T 


tih?Qk 2 } = ' 


] A\-j- «i{l} + 2r 2 ei{r 2 } - 87 £i{t 3 } + 8£i(t 4 }) + 

)Ar„ 2 (5r 2 ei{l} — 87€i{r} + 16 ci{t 2 }) + 6<r„ 4 £i{l} 


0 = T < ' 


£ 2 {l} - 9T>t 2 {r} + 15 r 2 e 2 {r 2 } - 12T£ 2 {t 3 } + 4« 2 {t 4 }) + 

_-4 2 cr„ 2 (7 T 2 c 2 { I } — 6 7<£ 2 {t} + 4£ 2 {r 2 }) + 6<r„ 4 £ 2 {l} -j- = t < T 


€{/* 2 4+l> = A\T* - 47 3 £{t} + 87 2 £{t 2 } - 87£{t 3 } + 4£{t 4 }) + A 2 cr n 2 (4T 2 - 87*0} + &0 2 }) + 4cr„ 4 


dQk 2 Qk+ 1> =< 


j^ 4 (-^— £l{l> + 2T 2 6\{t 2 } + 4£,{t 4 }) + 
I ^ 2 o-„ 2 (27' 2 £i{1} + 8£i{t 2 }) + 4cr„ 4 £]{l} 


0 S T < • 


A\ z5 j* £ 2 {1} - 207 3 £ 2 {t} + 26T 2 £ 2 {r 2 } - 167£ 2 {t 3 } + 4£ 2 {t 4 }) + 

^ 2 <T„ 2 (10r 2 £ 2 {l} - 16r£ 2 {r} + 8£ 2 {t 2 }) + 4<r„ 4 £ 2 {l} Z- § t < T 


e{Ik 2 C&+ 1> 


| — €i{l} — T^iO-} + 3T 2 ci{t 2 } — 47iei{r 3 } + 4€i{r 4 }) • 

] A 2 cr n 2 (3T 2 ei{l} — 4 T£\{t} + 8ei{r 2 }) -f 4o-* 4 £i{l} 


T < 


A*{^y~ £ 2 {l} - 97*£ 2 {t} + 157' 2 € 2 {t 2 } - 12r£ 2 {r 3 } + 4£ 2 {t 4 }) + 

^l 2 tr„ 2 (7r 2 £ 2 {l} — 127*20} + 8£ 2 {t 2 }) + 4<x„ 4 £ 2 {l} -j- = t < T 


t{ll+lQk 2 } = 



£]{1> - AiM + iT^d-r 2 } - 4J-£i{t 3 } + 4£,{t 4 » • 


A 2 <r„ 2 (5T 2 ei{l} — 27*10} + 4£i{t 2 }) + 6cr„ 4 £i{l} 


0 = t < -j- 


^ 4 (-^|— £ 2 {l} 22lL_ £2 { t } + _Z 1Z — e 2 {r 2 } - 157'£ 2 {t 3 } + 5£ 2 {t 4 }) + 

^W(13r 2 £ 2 {l} — 24r£ 2 {r} + 16£ 2 {t 2 }) + 6o-„ 4 £ 2 {l} § t < T 

\ *• 


(A- 14) 


(A-15) 
(A- 16) 


(A-17) 


(A-18) 


where, in the above equations 


dKr)} = j 

. T 

f(r)p(T)dT 

0 

(A-19) 

*lOO)} = 

rf 

I J f(~)p(r)dT 

(A-20) 

«2 OO)} = 

J T f(r)p(r)dr 

(A-21) 


where p(x) is the probabilty density function of the 
variable r. 


Appendix B 

5 g Derivation of the Mean and Variance of the AVOD 

Note that the calculations in this appendix incorpor- 
tate the results of Appendix C. The inphase and quadra- 
phase outputs are given, respectively, by Eqs. (13) and 
(14). The lock detector output X*= |Lj — |Qfc| . Let 
55 fik= e{ | lit | - | Q* | }• Then, 

Var(Xk) = e{l\Ik\ - \Qk\] 2 } - Pk 2 , (B-l) 

= till? + Qk 2 - 2 j 1 } - P-k 2 


Co«X k JC k+ j) = £{[1/*| - |fi*| - p.k][\Ik+j\ - (B-2) 

I Qk+j\ - Pk+fi 
= t{.\I k Ik+j[} + j QkQk~j \ } — 
t{\IkQk+j\} - t{\I k +jQk\} - PkPk+j 

The following equations were used to compute the 
variance of X k and the covariance of X* with X*+i: 
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\mi\AT- Ur + A7|} - e\{\Ur+ AT|}] 


0S '<f 


i[e 2 {\AT- Ur + AT|} - e 2 {\2AT - Ur + iV|}] -±- ^ t < T 


(B-3) 


t{\IkQk\} = 


<i 


mi\[AT+ Nt(k) + N 2 (k)][AT + N 2 (k) + N,(k + 1)]|} + 
t\{\\AT- 2 At +Ni(k) + N 2 (k)][-Ur + N 2 (k) + N x [k + 1)] | }] 

i[c 2 {\\AT + Ni(k) + N 2 (k)][AT + N 2 (k) + N x (k + 1)]|> + 
t 2 {\ [AT + Ni(k) + mk)][2AT -2 At + N 2 (k) + N\(k + 1)]|} + 

£ 2 {| [AT -2At + Ni(k) + N 2 (k)][UT - Ur + A r 2 (k) + N,(k + 1)]|} + 
€ 2 {| [AT -Ur + Ni(k) + N 2 (k)][^T + N 2 (k) + N\(k + 1)]|>] 


(B-4) 


o ^<7 


y s T < T 


«{|/*A+l|} = i[<?{\AT+N x (k) + JV 2 (i)|> + 

e{\[AT -Ur + N\{k) + N 2 (k)][AT - Ur + N x (k + 1) + N 2 (k + 1)) | } + 
2e{\ AT+ Ni(k) + N 2 (k)\}e{\AT - Ur + N X (k + 1) + N 2 (k + 1)|}] 




s' 

\[t\K\AT+ N\ + iV 2 |} + e\{\[Ur + Ni(k + 1) + AT 2 (*)] X 
[Ur + Ni(k + 2) + N 2 (k + 1)]|} + 

2e{\AT + N\ + JV 2 |}£i{|Z4r + N x +AT 2 |}] 

i[t 2 2 {\AT +N X + A 2 [} + e 2 {\[UT - Ur + AT, (* + 1) + N 2 (k)\ X 
[UT- Ur + N\{k + 2) + N 2 (k + 1)]|> + 

2«{|^r + AT] + N 2 \}t 2 {\UT - Ur + N x + iV 2 |» 


(B-5) 


(B-6) 


0§r<y 


fSr< T 


k[c 2 {\AT + Ni + Al 2 |} + €i{|t2^r + Ni(k + 1) + N 2 (k)] X 
[AT- Ur + N\(k + 2) + N 2 (k + 1)]|> + 

2e{\AT + N] + A-2|>6i{|2^t + N x + jV 2 |}] 

l/4[e 2 {\AT + N\ + AT 2 |} + £ 2 {| [AT- Ur + N\(k + 1) + N 2 (k)] X 
[UT -Ur + N\(k + 2) + N 2 (k + 1)]|} + 

^t{\AT + JVi + JV 2 |HMr - Ur + Ni + AT 2 |>] 


(B-7) 


0 — T < 2 


f^< r 


fi[ei{[[AT + N x (k + 1) + N 2 (k + 1 )][AT + N x (k + 1) + N 2 (k)]\} + 

*{\Ik+\Qk\} = < e l { tt Ar - 2Ar + N * k + D + N 2 (k + \)][AT+N x (k + 1) + W*||} + 

] t\i\[AT -Ur + N\(k + 1) + N 2 (k + \)][Ut + N\(k + 1) + JV 2 (/c)]|} + 
e\{\[AT + N x (k + 1) + N 2 (k + l)][Ur + N\(k + 1) + AT 2 (*)] ] >] 

-Ur + N\(k + 1) + N 2 (k + 1)][2^T - Ur + N x (k + 1) + 2V a (*)] | > + 
42 {|[^r + N\(k + 1) + N 2 {k + 1 )[[AT + N\(k + 1) + A^(*)] f >] 


(B-8) 


0Sr<y 


T = t < t 


where ei and e 2 are defined in Appendix A. The follow- 
ing functions were defined to obtain the results in Sub- 
section III.A of the main text. 

•fi(%)=«{ I (1 +cn\ + cn 2 )(l +cn 2 +CT3> [ > 

Fpl si =«{ | ( 1 + cn i + cn 2 )( cn 2 +C/13) | } 

Fjfyi | (1 — 2a + cn t + cn 2 )( 1 — 2a + <773 + c«4> | } 

GiCtIs) = ei{ i (1 — 2a + cni + cn 2 )(2u + cn 2 + cn 2 ) | } 

G 2 (T)j)=ei{ I (1 -2u+cn\ +cn 2 )( 1 +cn 2 +cn 2 ) \ } 


G 3 (%) =e 1 { | (2u + cn 1 + crt 2 )(2u + cn 2 + cni) I ) 
c?4(%) =£ 1 { | (2a + cn 1 + cn 2 )( 1 - 2u + cn 2 + cni) I } 

60 -ffl(%) =« 2 { I (1 +cn\ + cn 2 )( 2 -2u+cn 2 +cni)\} 

H 2 (t]s) =«{ [ (1 — 2u + cn 1 + cn 2 )(2 — 2u + cn 2 +cn 2 ) | > 

^3(7 lx) =e{ | ( 1 — 2u + cn 1 + cn 2 )( 1 + cn 2 + cn 2 ) \ } 

# 4 ( 1 ) 1 ) =e{ | (2 — 2 a + cn \ + cn 2 ){2 — 2 a + cn 2 + cni) I ) 

# 5 ( 1 ) 1 ) =£{ I (1 - u + cn 1 + cn 2 )(2 - 2a + m 3 + cni) I } (b-9) 
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where c = 1 /( 2 %/%), the n,’s are normal independent 
random variables with zero mean and unit variance, and 
u is a uniform random variable in the range [0, 1]. These 
functions have been computed as follows: In the F func- 
tions, the expectation with respect to u is carried over 5 
the entire region [0, 1], while in G and H functions, the 
expectation is carried over [ 0 , £] and [i, 1 ], respectively. 


t\{\AT- IAt + iVi + N 2 II = 



Vs ) + 


(C-6) 


Appendix C 

Derivation of the Mean for the Random Variable 10 
|n+bT+c| 



+ 



exp(-%) 


Let n be a normal random variable with zero mean 
and variance cr 1 , t a uniform random variable over (0, 
T), and c any constant. Then, 


«{|n + c|} = 


N 2tt(T“ 



\n + c|exp 



(C-l) 


Also, by simple manipulation, it can be shown that 


15 (C-7) 

«l{| AT- 2At + Ni + N 2 \} = <l 2 {\AT-2At + Ni + AT 2 |} 

= 4- N\ 4- 1V 2 |} 

= e 2 {|2 AT - 2At + Ni + 

N 2 \} 


The above integral can be easily evaluated by breaking 
it into two integrals over the two regions (— 00 , — c) 
and (— c, co) to obtain 


20 


e{ l' ! + ei} = \|^ exp (-i') +cer/ (Y^) 


Appendix D 

Derivation of the Mean and Variance of the SQNOD 

The samples Land Q*are given by Eqs. (33) and (34). 
(C-2) 25 The moments of X*, are given by Eqs. (A-l) through 
(A-3). Using Eqs. (D-7) through (D-10) with r= 0 and 
Eqs. (A-2) and (A-3) yields the in-lock variance. 
Namely, 


For a fixed t, one can write 
c{[n + br + c|} = 


30 


(C-3) 


VartJC, t) = + -y A 2 T 2 cr„ 2 + 4c r„ 2 


P-1) 


|^ esp (_l^) +( . + c>r 0) 

Unconditioning over r yields 

(C-4) 


The covariance can be shown to be, for ji£ i ( 

35 cov(XbXk+j)=0 P-2) 

Similarly, with r modelled as uniform over [0, T/2] in 
Eqs. (D-7) through (D-10), the out-of-lock moments are 
found. Hence, 


£l{|« + br + c|} = 4 t 
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VarJiXk) = + y A 2 T 2 a„ 2 + 4<r„ 4 


P-3) 


(br + c)erf 


which, after integration by parts, leads to 


(w) 


The out-of-lock covariance can be shown to be, for 
* 45 jSl, 


CovJiX k JC k+ j) = 


P-4) 



By applying the above expression, one gets 


The following equations are used in computing the 
variance of X*: 
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, a 2 t 2 


<K/* 2 } = 


+ <r„ 2 )ei{l> 


0 S r < 


tiQk 2 ) = 


A 2 (— —— «2(l} - 2T<L2{t} + 2£2{t 2 >) + cr„ 2 e 2 O} %t <—■ 

f 

A 2 (-~— «,{!} + 2«i{t 2 }) + <r n 2 €i{l} 0 S r < -j- 


T 2 
8 

A 2 ! 2 


+ Or„ 2 )« 2 {l} 


£-<t 


(D-5) 


(P-6) 


+ T^ 2rV » 2 + 3 ^ 4 ) e i« 

e{/* 4 } = | A\- 

JA . 

4— r 2 € 2 {0 - 12 T€ 2 {t} + 12€ 2 {t 2 }) + 3<r n 4 e 2 {l} 


16 ~ r 2 

«2-Cl} - 47 6 £2 {t) 4- 12r 2 « 2 {T 2 } - 16762 -Ct 3 } + 8C2{t 4 }) 

A 2 cr„ 2 (- 15 


»Sr<f 


T_ 

4 


i T < -=- 


£{& 4 } = 


A 4 (-^~ £l(l} + 8€,{r 4 }) + A 2 cr 2 (- r 2 ti{l} + !26 i{t 2 }) + 3tr„ 4 ei{l} 0 ^ r < ~ 


. ^T 4 
■ 16 


- A 2 T 2 <x„ 2 + 3<r„ 4 )t2{l} 


T_ 

4 


(D-7) 


(D-8) 


Appendix E 

Derivation of the Mean and Variance of the AVNOD 

By following the same procedure as in Appendix B 
and by using 


u is a uniform random variable in the range (0, 1). In Z, 
the expectation with respect to u is over the range [0, $]. 

Appendix F 

Derivation of the Mean and Variance of the SPED 




h{|4f- + JV|} 


0 = T < ' 


i[£2{| 4 y : - + A|} + £2{|-47'-2^T + Ar|}] |-£r<j 


(E-l) 


£{|fit|> = 


I H£li|- 4 Y- + A r |} + filler + 1V|}] OS t<£- 
U>i\^- + N\} 


(E-2) 


and 


s' 

~ €j{l} + 2A 2 €\{t 2 } + 2cr„ 2 €i{l} - 

J 5/4ei 2 {|-4p- + 2V|} - }£1 2 {|2^It + tV|> - 

Var(\Ik\ - \Qk\) = < , r , AT , 

| I«{| — — + AT| Jei-CIZ^-r + JV[} 


(E-3) 


o^<{ 


T ~ - «2{1) - 2A 2 Tt2{r) + 14 2 £2{t 2 } + 2cr„ 2 e2{l} - 
5/4£2 2 f| 4 Y- + Nl}-i£2 Z {|-4T-MT + lV|>- 

+ N\}i 2 {\AT - 2At + N\} S T < y- 


one obtains Eqs. (40) through (45) after using the results 
of Appendix C and lengthy manipulations. The function 
Z in Eq. (45) is defined as 

Z(rts)^;€\{\(u+cni)(u+cn 2 )\} (E-4) 65 

The samples I*, and Qk are given by Eqs. (47) and 
where c = 1/(2Vt] x ), the n,’s are normal independent (48). The output of the lock detector is Xk=lkQk- It is 
random variables with zero mean and unit variance, and straightforward to show that 
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•Utf- A f - A \ T <r) 

(F-l) 

. T 2 t{r 2 } 

V CM) 5 

; 6 4 + 2 

) 

<r„ 4 + A 2 o-„ 2 

- Te{r} + 2«{r 2 } 'j ^ 


(F-3) 

Co Wk+ j)j 4 

y-^and/i; 1 


0 T = 0 15 


Equations (49) and (50) follow after letting r=0 in Eqs. 
(F-l) through (F-3). Equations (51) and (52) follow by 
letting r be uniform over [0, T/2]. 

We claim: 

1. A symbol lock detector for use in coherent digital 
communication systems after synchronization of a car- 
rier and any subcarrier has been achieved, where sym- 
bol data, dk, takes on binary values with equal probabil- 
ity, and symbol interval, T, is known, comprising 

a single integrator for continuously integrating an 
incoming signal to form in sequence pairs of nono- 
verlapping first and second integrals over succes- 
sive intervals, each integral selected to be over an 
interval T/2, 

means for delayinq said first inteqral of each of said 
pairs of integrals for said interval, T/2, in transmit- 
ting said first inteqral of each of said pairs of inte- 
grals, 

means for multiplying said first integral delayed by 
delay means and said second integral of each of 
said pairs of integrals to produce a value Xfo 

means for multiplexing said first integral of each of 
said pairs of integrals to said delay means and for 
multiplexing said second integral of each of said 
pairs of integrals directly to said means for multi- 
plying, thereby providing for multiplication of said 
first integral of each of said pairs of integrals by 
said second integral of each of said pairs of inte- 
grals, 

means for averaging said value X* over a number of 
symbol intervals to produce a statistic value y, and 

means for thresholding said statistic value y at a level 
e to produce a symbol lock decision when said 
statistic value y exceeds said level e. 

2. A symbol lock detector for use in coherent digital 
communication systems after synchronization of a car- 
rier and any subcarrier has been achieved, where sym- 
bol data, dfo takes on binary values with equal probabil- 
ity, and symbol interval, T, is known, comprising 

a single integrator for continuously integrating an 
incoming signal to form in sequence pairs of nono- 
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verlapping first and second integrals over succes- 
sive intervals, each integral selected to be over an 
interval T/2, 

a single means for squaring each of said first and 
second integrals of each of said pairs of integrals in 
sequence, thereby producing first and second 
squared integrals, 

means for delaying said first squared inteqral with 
respect to said second squared inteqral by said 
interval T/2, 

means for summing said first squared integral delayed 
by delay means and said second squared integral to 
produce a value X*, 

means for multiplexing said first squared integral to 
said delay means and said second squared integral 
to summing means, thereby forming a sum of said 
first squared integral delayed by said delay means 
and said second squared integral that is said value 
X fe 

means for averaging said value X* over a number of 
symbol intervals to produce a statistic value y, and 
means for thresholding said statistic value y at a level 
£ to produce a symbol lock decision when said 
statistic value y exceeds said level e. 

3. A symbol lock detector for use in coherent digital 
communication systems after synchronization of a car- 
rier and any subcarrier has been achieved, where sym- 
bol data, dk, takes on binary values with equal probabil- 
ity, and symbol interval, T, is known, comprising 
a single integrator for continuously integrating an 
incoming signal to form in sequence pairs of nono- 
verlapping first and second integrals over succes- 
sive intervals, each integral selected to be over an 
interval T/2, 

a single means for forming the absolute value of each 
of said first and second integrals of each of said 
pairs of integrals in sequence, thereby producing 
first and second absolute-value integrals, and 
means for delaying said first absolute-value inteqral 
with respect to said second absolute-value inteqral 
by said interval T/2, 

means for summing said first absolute-value integral 
delayed by delay means and said second absolute- 
value inteqral to produce a value X/t, 
means for multiplexing said first absolute-value inte- 
gral to said delay means and said second absolute- 
value integral directly to said summing means, 
thereby forming a sum of said first absolute-value 
integral delayed by said delay means and said sec- 
ond absolute-value integral that is said value X*, 
means for averaging said value X* over a number of 
symbol intervals to produce a statistic value y, and 
means for thresholding said statistic value y at a level 
S to produce a symbol lock decision when said 
statistic value y exceeds said level 8. 

* * * * * 
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